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EQUILIBRIUM-VAPOR  CELL  FOR  QUANTITATIVE  INFRARED 
ABSORBANCE  MEASUREMENTS 


1 .  INTRODUCTION 

The  detection  and  quantification  of  volatile  organic 
compounds  (VOCs)  at  hazardous  waste  sites  are  important  in 
determining  the  necessary  steps  for  remediation.  Open-path 
Fourier  transform  infrared  (FTIR)  spectrometry  provides  an 
effective  method  of  monitoring  for  these  VOCs.‘’^  A  variety  of 
approaches  for  detection  of  VOCs  with  open-path  FTIR  spectrometry 
are  documented.^”* 

This  investigation  furnishes  a  vapor-cell  method  for 
generation  of  known  vapor  concentrations.  This  vapor-cell 
approach  differs  from  the  traditional  methods  of  vapor  generation 
for  static  gas  cell  measurements.’’*  The  equilibrium-vapor  cell 
(EVC)  relies  on  the  maintenance  of  a  constant  vapor  concentration 
by  establishment  of  equilibrium  with  the  continuous  circulation 
of  an  air-vapor  mixture  through  a  fixed  composition  aqueous 
solution.  Aqueous  solutions  are  used  for  convenience.  These 
solutions  also  permit  a  reliable  assessment  of  water  vapor 
concentration,  as  well  as  VOC  solute  vapor  pressure.  Aqueous 
solutions  necessitate  that  the  VOC  be  completely  miscible. 

Binary  aqueous  solutions  of  methanol,  ethanol,  isopropanol, 
acetone,  and  methyl  ethyl  ketone  (MEK)  are  used. 

Experimental  determination  of  the  vapor  concentrations 
of  non-ideal  aqueous  solutions  is  extremely  difficult.  The  EVC 
method  relies  on  the  accurate  weighing  in  solution  preparation 
and  on  the  proper  temperature  equilibration  in  using  the  prepared 
aqueous  solution.  Inputs  of  liquid  composition  and  solution 
temperature  permit  determination  of  the  vapor  concentration  for 
the  solute  and  solvent  with  the  Wilson  equation.®  The  Wilson 
equation  is  accurate  to  within  2%  for  the  systems  studied.  This 
requires  a  realistic  assessment  of  the  Wilson  coefficients  from 
the  available  literature  vapor  data.*®  The  EVC  absorbance 
measurements  are  made  with  a  conventional  MIDAC  (Irvine,  CA)  FTIR 
spectrometer.  The  absorbance  measurements  are  shown  adherence  to 
Beer's  Law  over  a  range  of  vapor  pressures.  The  slope  of  the 
absorbance  versus  the  vapor  concentration  is  the  absorptivity- 
pathlength  product.  The  absorptivity  is  obtained  by  dividing 
this  product  with  the  cell  pathlength.  The  experimentally 
determined  absorptivity  is  compared  with  those  from  other 
laboratories  for  five  solute  vapors. 


2.  EXPERIMENTAL  METHODS 

The  FTIR  spectrometer  is  equipped  with  a  liquid 
nitrogen  cooled  Hg:Cd:Te  detector  (model  no.  HCT-21A, 
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serial  no.  P-14319,  Infrared  Associates,  Incorporated,  Cranbury, 
NJ) .  The  laboratory  spectrometer  is  a  MIDAC,  model  101280 
(serial  no.  106,  MIDAC  Corporation,  Irvine,  CA) .  The  laboratory 
spectrometer  detector  has  a  spectral  range  of  400  to  4000  cm'*. 

Interferogram  data  was  collected  with  Spectra  Calc  at 
2  cm'*  spectral  resolution  for  all  vapors  generated  with  aqueous 
solutions,  as  well  as  0.56  cm*  resolution  for  only  the  methanol 
vapor.**  Interferograms  of  100  coadditions  were  collected  with 
the  laboratory  spectrometer.  Analysis  of  the  interferogram  and 
spectral  data  was  performed  in  Spectra  Calc. 

The  equilibrium-vapor  cell  consists  of  readily 
available  components.  Figure  1  is  a  schematic  of  the 
equilibrium-vapor  cell.  The  circulating  pump  model  81-012 
(Universal  Electric  Company,  Owosso,  MI) ,  demountable  10-cm  gas 
cell,  model  7200  (Buck  Scientific,  East  Norwalk,  CT)  with  zinc 
selenide  windows,  and  bubbler  are  interconnected  with  tygon 
tubing.  The  bubbler,  which  contains  125  mL  aqueous  solution,  is 
immersed  in  a  constant  temperature  bath  model  RTE-111  (Neslab 
Instruments  Incorporated,  Newington,  NH) .  The  recycling  of  the 
air-vapor  mixture  through  the  aqueous  solution  achieves 
equilibrium  in  the  gas  cell  after  approximately  10  min.  The 
bath  is  set  at  a  temperature  of  22  «C,  which  is  below  the  ambient 
room  temperature.  This  avoids  condensation  of  vapors  in  the 
connection  tubes,  pump,  or  demountable  cell.  The  bath  temper¬ 
ature  is  monitored  with  a  calibrated  negative  temperature 
coefficient  thermistor. 

Aqueous  solutions  were  prepared  by  weight  percent. 
Weights  of  solution  components  were  made  on  a  model  AE200 
top  loading  balance  (Mettler  Instruments  Corporation,  Hightstown, 
NJ) .  Distilled/deionized  water  was  used  in  solution  preparation. 
The  sequence  of  solutions  used  in  the  absorbance  measurements 
proceeded  from  most  dilute  to  most  concentrated.  This  avoided 
any  potential  problems  with  vapor  absorption/desorption  in  the 
equilibrium- vapor  cell  apparatus. 


3.  RESULTS  AND  DISCUSSION 

The  dependence  of  the  IR  absorbance  on  the  partial 
vapor  pressure  of  five  organic  solvents  is  examined  for  methanol, 
ethanol,  isopropanol,  acetone,  and  MEK.  The  partial  vapor 
pressures  for  all  but  three  measurements  range  from  5  to 
50  millibars  at  22  ®C.  Various  concentrations  of  binary  aqueous 
solutions  provide  seven  evenly  distributed  vapor  pressures.  The 
liquid  mole  fractions  and  vapor  pressures  at  22  ®C  of  the  binary 
solutions  are  tabulated  in  Table  1. 

Accurate  determination  of  the  vapor  pressure  is 
essential  to  the  quantitative  evaluation  of  the  IR  absorbances. 
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All  of  the  binary  solutions  considered  in  this  study  exhibit  non¬ 
ideal  solution  behavior.  Therefore,  the  empirical  evaluation  of 
the  activity  coefficients  of  the  solution  components  becomes 
necessary  to  account  for  the  excess  solution  free  energy  of 
mixing.  Rather  than  rely  on  a  single  data  set,  the  approach  of 
this  investigation  follows  a  modified  version  of  Hirata  et  al.'® 
to  obtain  the  best  pair  of  Wilson  coefficients  for  each  binary 
solution.  Based  on  the  analysis  of  the  various  empirical  methods 
for  evaluation  of  the  activity  coefficients,  Hirata  et  al. 
maintained  that  Wilson's  equation  is  the  most  reliable  for  a 
broad  range  of  systems  including  those  in  this  study. 

The  Wilson  equation  accounts  for  the  excess  solution 
free  energy  of  mixing  by  extending  the  Flory  and  Huggin's 
theoretical  equation  for  athermal  solutions.®  The  relative 
probabilities  of  dissimilar  neighboring  molecules,  within  a  local 
volume  fraction,  permits  derivation  of  the  Wilson  coefficients. 
These  coefficients,  along  with  solution  mole  fractions,  allow  a 
determination  of  the  excess  free  energy  of  mixing.  The  Wilson 
equation  for  determining  the  activity  coefficient  is  subsequently 
derived  from  the  excess  free  energy  of  solution  mixing.  The 
Wilson  equation  for  the  ith  component  in  a  binary  system  (i,j)  is 
described  by  equation  1 . 


Xj^Wj.X, 


(1) 


X  is  the  mole  fraction  in  solution.  The  Wjj  and  Wjj  are  the  Wilson 
coefficients.  An  advantage  of  this  equation  is  that  the  general 
equation  for  three  (i,j,k)  or  more  components  requires  only  the 
corresponding  pairwise  coefficients  (i.e.,  Wii.,Wkj,Wki,Wkj)  , 
eliminating  the  need  for  triplet  coefficients. 

Once  the  system  Wilson  coefficients  are  known,  the 
vapor  pressure,  Pj,  of  the  individual  components  are  computed 
from  the  activity  coefficients,  ji,  as  shown  in  equation  2. 

Pi=yi^iPi  (2) 


P°  is  the  vapor  pressure  of  the  pure  component  at  the  solution 
temperature.  The  empirical  Antoine  equation  provides  the  pure 
vapor  pressures  in  equation  3. 


In  pj=  (In  10) 


(3) 
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Aj,  Bj,  and  C;  are  the  experimentally  determined  Antoine 
coefficients  for  the  ith  solution  component,  and  t  is  the 
solution  temperature  in  Celsius.  For  binary  solutions,  1 
represents  the  organic  solute  and  2  represents  the  solvent  water. 

The  original  literature  data  cited  by  Hirata  et  al. 
established  the  pairs  of  Wilson  coefficients  for  the  five  binary 
systems.  These  data  files  range  from  as  few  as  six  references 
for  water-MEK  (50  data  points)  to  as  many  as  22  references  for 
water-ethanol  (298  data  points) .  The  original  vapor  pressure 
data  for  each  of  these  files  are  either  isobar ic  or  isothermal. 
Isobaric  indicates  that  the  vapor  mole  fraction  and  temperature 
are  in  equilibrium  with  the  solution  mole  fraction  at  a  specified 
constant  total  pressure.  Isothermal  refers  to  the  condition  of 
the  vapor  mole  fraction  and  total  pressure  in  equilibrium  with 
the  solution  mole  fraction  at  a  specified  constant  temperature. 
Besides  the  binary  solution  vapor  pressures,  the  vapor  pressures 
of  the  pure  components  are  also  needed  in  computing  the  Wilson 
coefficients.  Following  the  approach  of  Hirata  et  al.,  these 
coefficients  are  evaluated  using  two  distinct  sets  of  Antoine 
coefficients  for  all  but  the  MEK  component.  Table  2  summarizes 
the  sets  of  Antoine  coefficients  for  the  systems  considered. 

There  is  not  an  apparent  correlation  or  consistency  in  which 
Antoine  coefficient  set  is  applied  by  Hirata  et  al.  for  a  given 
binary  solution.  In  a  few  instances,  selecting  one  or  the 
alternate  set  of  Antoine  coefficients  yields  unrealistic  (i.e., 
negative)  Wilson  coefficients  that  are  unacceptable.  The 
differences  in  the  vapor  pressures  of  the  pure  components 
predicted  by  each  Antoine  coefficient  set  are  small  and  appear  to 
furnish  a  better  vapor  pressure  fit  for  either  the  high  or  low 
ends  of  the  liquid  temperature  range.  Rather  than  an  arbitrary 
selection  of  Antoine  coefficients,  the  present  study  recalculates 
the  pairs  of  Wilson  coefficients  for  each  literature  citation 
(data  file) ,  using  all  four  combinations  of  Antoine  coefficients 
(i.e.,  water-organic,  W-0:  1-1,  1-2,  2-1,  and  2-2).  The  pair  of 

Wilson  coefficients  yielding  the  best  fit  is  determined.  The  set 
of  Antoine  coefficients  that  resulted  in  the  best  fit  Wilson 
coefficients  are  listed  in  the  last  two  columns  of  Table  3.  The 
criteria  for  best  fit  of  the  Wilson  coefficients  is  the  smallest 
average  deviation  of  the  vapor  mole  fraction  given  in  equation  4 . 


6y  = 


Y^{calc)  -  Y^(obs)\ 


(4) 


The  average  deviation  for  the  temperature,  ST,  is  similarly 
defined.  These  figures  of  merit  are  consistent  with  those  used 
by  Hirata  et  al.  and  allow  an  additional  verification  of  the 
vapor  pressure  computations.  The  percent  variation  in  the 
average  deviation  of  the  vapor  mole  fraction,  %(JY,  for  the  water- 
methanol  system  (i.e.,  100*(SY)  is  tabulated  in  column  6  of 
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Table  3;  whereas  the  average  deviation  for  the  temperature,  ST, 
is  listed  in  column  7. 

The  computation  of  the  pair  of  Wilson  coefficients  for 
a  given  data  file  minimizes  the  residual  sum  of  squares  for  the 
excess  Gibbs  free  energy  of  solution: 

SQ  =  J^z\,  (5) 


where 


•^i  =  E  -^i  In  Yi(oi)s)  -  E  InYi(caJc) , 


(6) 


with 


Yi(obs) 


Uij 


(V) 


and  In  7j(calc)  is  calculated  from  the  Wilson  equation  (i.e., 
equation  1)  .  Yj  is  the  vapor  mole  fraction  and  Pj  is  the  total 
vapor  pressure  in  equation  7.  Substitution  of  equation  1  and 
equation  7  into  equation  6  produces  a  non-linear  expression  for  a 
two  component  system,  namely  equation  8; 


= 


/  y  \  /  y  \ 

\-^l)  \  )  \  "*2  / 


(8) 


The  minimization  of  SQ  is  performed  by  setting  (3SQ/5Wi2)=0  and 
(5SQ/W2i)=0.  This  gives  two  difference  equations  in  AW12  and  AW21. 

A  f^l2  E  ^  ^2iE  iP^i  =E  (  9 


A  E  ^  ^  ^2iE  =E 
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where 


Dl,= 


(11) 


D2,= 


\x,^x^^x,]. 


(12) 


Equations  9  and  10  that  are  defined  by  equations  11  and  12  are 
solved  simultaneously  for  AWjj  and  AWji  by  the  proper  selection  of 
values  for  W12  and  W21.  Initial  values  of  0.5  are  selected  for  the 
Wilson  coefficients.  The  computation  of  Wilson  coefficients  is 
subsequently  iterated  by  setting  Wij(new)=Wij(old) -AWjj.  This 
iteration  process  continues  until  the  differences  are  below  a 
desired  threshold  value  (e.g.,  | AWjj |  + 1 AW21 1 <10-’  to  be  consistent 

with  Hirata's  computational  approach). 

The  calculation  of  the  Wilson  coefficients  for  all  data 
files  of  a  specific  two  component  systems  produces  one  pair  of 
Wilson  coefficients  for  each  file.  Hirata  et  al.  have  observed 
that  variation  of  W12  with  respect  to  W21  for  a  given  data  file 
results  in  a  minimum  in  SQ  that  shifts  with  the  value  of  W21. 

Thus,  a  single  minimum  exists  along  this  hyperbolic  path. 

Although  Hirata  et  al.  do  not  address  the  selection  of  the  best 
single  Wilson  coefficient  pair  to  represent  all  the  data  files 
given  for  a  two  component  system,  this  investigation  demonstrates 
all  the  systems  possess  a  similar  single  minimum  for  SY .  The 
pair  of  Wilson  coefficients  for  this  minimum  are  close  but  not 
identical  to  the  mean  values  calculated  over  the  array  of  pairs 
computed  from  the  data  files.  Table  3  shows  these  results  in 
columns  4  and  5.  The  actual  ”best  pair"  of  Wilson  coefficients 
are  found  by  a  grid  search  spanning  the  entire  array  of  pairs. 
Figure  2a  presents  a  global  %(5Y  error  graph  for  W21  versus  W2, , 
whereas  Figure  2b  is  the  projected  two  dimensional  contour  plot 
of  the  minimum  hSY  errors  derived  from  the  grid  search.  The  hSY 
value  that  is  designated  by  the  dash  lines  in  Figure  2b  locates 
the  minimum  for  the  aqueous  methanol  system.  The  Wilson 
coefficient  values  at  the  minimum  %(SY  are  given  in  Table  3  as 
the  best  values.  Table  3  provides  a  summary  of  the  computed 
results  for  the  water-methanol  system.  Columns  8,  9,  and  10 
(Table  3)  list  the  percentage  for  the  individual,  mean,  and  best 
estimate  of  the  standard  error  in  the  vapor  mole  fraction, 
respectively.  The  Wilson  coefficients  for  the  water-methanol 
system  are  reported  in  the  last  two  rows  of  Table  3,  which  are 
labeled  MEAN  and  BEST.  In  calculation  of  the  mean  and  best 
Wilson  coefficients,  the  W12  coefficient  associated  with  reference 
file  number  523  is  rejected  as  an  outlier  by  the  Grubb's  test'^'^ 
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at  the  99.5%  confidence  level.  Thus,  both  Wilson  coefficients 
for  reference  file  number  523  must  be  rejected.  Table  4  supplies 
a  synopsis  of  the  results  for  the  five  binary  systems  studied 
with  the  mean  and  best  Wilson  coefficient  values. 

The  pair  of  Wilson  coefficients  established  for  each 
system  allow  calculation  of  the  vapor  pressures  at  22  ®C  for  the 
seven  solutions  of  each  binary  system,  based  on  the  solution  mole 
fractions  with  the  exception  of  only  four  concentrations  for  the 
MEK -water  system.  The  mole  fractions  of  each  solution  are 
computed  directly  from  the  measured  weights  of  the  components. 
These  weights  have  been  determined  to  be  within  0.01  g  for  total 
solution  weights  of  around  170  g.  Table  1  summarizes  the 
compositions  and  vapor  pressures  that  are  used  in  accessing  the 
adherence  of  the  associated  IR  absorbances  to  Beer's  Law. 

Infrared  interferograms  are  acquired  with  the 
equilibrium  vapor  cell  for  32  aqueous  solutions  equilibrated  at 
22  °C.  The  solution  compositions  are  listed  in  Table  1  by 
percent  solute  mole  fraction  in  the  solvent  and  the  associated 
solute  vapor  pressure  developed  over  the  solution.  The  methanol 
vapor  provides  the  opportunity  to  investigate  the  onset  of 
rotational  fine  structure  at  1  cm'*  spectral  resolution  for  the 
P  and  R  branches  near  1014  and  1054  cm'*,  respectively.  The  PQR 
branch  contour  centered  at  103  3  cm'*  is  devoid  of  other  spectral 
features  in  the  800  to  1200  cm*  region.  Integration  of  the  PQR 
absorbance  band  contour  and  use  of  the  path  averaged  concen¬ 
tration  permits  a  comparison  of  the  absorbance  measurements 
for  this  investigation  to  those  of  other  laboratories.  The 
absorbance  measurements  for  methanol  are  made  at  0.56,  1,  2, 
and  4  cm*  spectral  resolutions.  These  measurements  allow  an 
evaluation  of  the  interferogram  apodization  effect  on  the 
absorbance  measurements.  The  remaining  solute  absorbances  are 
determined  from  interferograms  with  a  2  cm'*  spectral  resolution. 
Select  absorbance  bands  in  the  spectral  region  from  800  to 
1400  cm*  are  used  in  calculation  of  the  absorptivities .  The 
absorptivities  that  are  obtained  with  the  EVC  approach  are 
compared  to  those  of  other  investigators. 

Methanol  vapor  pressures  are  varied  from  5.7  to 
37.9  torr,  which  corresponds  to  a  path  averaged  concentration 
from  758  to  5039  ppm-m  with  a  10-cm  pathlength  cell.  The  corre¬ 
sponding  water  vapor  pressures  range  from  19.1  to  15.4  torr  with 
the  path  averaged  concentrations  ranging  from  2540  to  2048  ppm-m. 
The  methanol  absorbances  for  a  methanol  vapor  pressure  are  about 
1000  times  stronger  than  those  of  water  at  2000  ppm-m  in  the  800 
to  1200  cm*  spectral  region.*^  The  blank  cell  serves  to  ratio  out 
this  contribution  from  the  water  vapor  absorbances.  The 
absorbance  measurements  for  the  PQR  bands  at  0.56  cm'*  spectral 
resolution,  using  a  box  car  apodization,  are  plotted  as  a 
function  of  vapor  pressure  in  Figure  3.  The  open  symbols 
represent  the  least  squares  calculated  absorbances  and  the  solid 
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symbols  denote  the  measured  values.  The  P  and  R  bands  adhere 
well  to  Beer's  Law  with  absorbance  versus  vapor  pressure  slopes 
of  0.03200  ±  0.00029  and  0.04034  ±  0.00022,  respectively.  On  the 
other  hand,  the  Q  branch  shows  a  marked  deviation  from  Beer's  Law 
with  a  slope  of  0.0530  ±  0.0017.  This  deviation  from  Beer's  Law 
occurs,  because  the  spectrometer  resolution  of  0.56  cm'*  is 
inadequate  to  resolve  the  narrow  Q  band. 

The  Beer's  Law  slopes  for  the  methanol  P  and  R  branches 
are  shown  as  a  function  of  spectral  resolution  and  apodization  in 
Table  5.  The  lower  spectral  resolutions  are  generated  by  the 
appropriate  truncation  of  the  0.56  cm'*  resolution  interferograms . 
The  effect  of  apodization  does  not  generally  become  pronounced 
until  a  spectral  resolution  of  1  cm'*  or  greater  is  reached.  This 
has  also  been  confirmed  for  all  other  solute  vapor  spectra 
collected  at  2  cm'*  spectral  resolution  in  this  study.  At  1  cm'* 
spectral  resolution,  the  magnitude  of  absorbance  versus  vapor 
pressure  slopes  follow  the  expected  trend  in  the  apodized 
absorbance  spectra  of  box  car  >  medium  Norton  Beer  >  triangular. 

The  integrated  absorbance  area  of  a  spectral  band  is 
independent  of  the  spectrometer's  line  shape  function.*^  There¬ 
fore,  the  absorbance  areas  permit  a  comparison  of  experimental 
results  between  spectrometers  that  are  operated  at  different 
spectral  resolutions.  A  comparison  of  methanol  absorbance 
spectra  from  this  investigation  to  those  of  three  other 
laboratories  is  made  with  Figure  Figure  4  shows  a  plot 

of  PQR  band  absorbance  area  versus  path  averaged  concentration 
(e.g. ,  ppm-m) .  The  limits  of  band  integration  are  established  by 
computing  the  best  estimate  of  the  standard  deviation  (BESD)  over 
two  180  cm'*  intervals  positioned  on  each  side  of  the  PQR  band 
contour.  The  lower  interval  begins  at  700  cm*  and  the  upper 
starts  at  112  0  cm'*.  The  BESDs  for  these  intervals  are  calculated 
with  successive  differences  as  shown  in  equation  13. 


BESD  = 


n-l 


i^l 


2[n-l] 


(13) 


The  Pi  denotes  a  specified  wavenumber,  and  A(j'i)  represents  the 
absorbance  at  that  specified  wavenumber.  The  application  of  the 
successive  difference  method  for  calculation  of  the  BESD  for 
spectra  is  documented  by  Mark  and  Workman.**  The  first  absorbance 
value  approaching  the  PQR  band  that  exceeds  three  times  the 
calculated  BESD  is  designated  as  an  integration  limit.  There¬ 
fore,  the  integration  limits  are  effectively  defined  in  terms  of 
the  absorbance  spectra  noise.  Integration  of  the  PQR  absorbance 
band  is  subsequently  performed  with  the  trapezoidal  integration 
routine  available  in  Spectra  Calc.  The  integration  results  for 
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the  methanol  PQR  band  are  shown  in  Figure  4.  The  integrated 
absorbances  demonstrate  reasonable  agreement  with  the  exception 
of  the  measurement  in  reference  16. 

A  comparison  of  absorbances,  at  the  same  spectral 
resolution,  is  required  to  compare  the  absorptivities  at  a 
particular  wavenumber.  The  absorbance  spectra  from  references 
14,  17,  and  16  possess  resolutions  of  0.125,  0.25,  and  2  cm'*, 
respectively.  Therefore,  the  absorbance  spectra  for  references 
14  and  17  are  adjusted  to  a  2  cm'*  resolution  for  purposes  of 
comparison  by  introduction  of  fractional  values,  p  (i.e.,  either 
0.0625  or  0.125)  into  equation  14. 


filter  #s, 0, l,p, l,p, 0, 1, 0 


(14) 


The  expression  given  in  equation  14  is  the  Spectra  Calc  Fourier 
filtering  function  that  reduces  the  spectral  resolution  of  the 
absorbance  source  spectrum  by  1/p. 

The  absorbance  spectra  at  2  cm*  resolution  are  used  in 
calculation  of  the  absorptivities.  If  only  a  single  measurement 
is  available,  then  Beer's  Law  is  solved  algebraically  for  the 
absorptivity,  a.  When  three  or  more  measurements  are  available, 
the  absorptivity  is  calculated  from  the  least  squares  slope  of 
absorbance  versus  concentration.  The  gas  concentration  in  grams 
per  cubic  meter  is  determined  from  the  vapor  pressure  in  torr  by 
the  ideal  gas  law  that  is  shown  in  equation  15. 

J  g\_  16.034MP,,^  (15, 

UV  T 


T  is  the  vapor  temperature  in  degrees  Kelvin;  M  is  the  molecular 
weight  of  the  vapor  solute  in  grams/mole /  is  the  vapor 
pressure  in  torr;  and  16.034  is  the  reciprocal  of  the  ideal  gas 
law  constant  in  units  of  mole  Kelvin/cubic  meter  torr  (mol*K)/ 
(m^*torr) .  The  final  units  for  the  absorptivity  are  square 
meters/milligram  (m^/mg) .  The  absorptivities,  a,  are  calculated 
for  a  2  cm*  spectral  resolution  and  are  summarized  for  all  solute 
vapors  studied  in  Table  6. 

The  resultant  absorptivities  in  Table  6  show  signi¬ 
ficant  variation  between  the  laboratories.  The  experimental 
precision  in  the  EVC  measurements  for  the  bands  selected  in 
Table  6  is  as  follows:  Methanol/P,  ±0.8%,  /Q,  ±4%,  and  /R,  ±0.5%; 
Ethanol/  ±2%;  Isopropanol/  ±2%;  Acetone/  ±0.5%;  and  MEK/  ±2%. 

Only  the  1072  cm*  band  of  isopropanol  for  reference  14  falls 
within  the  EVC  experimental  precision.  The  remaining  isopropanol 
absorptivities  for  reference  14  are  within  about  7%  of  the  EVC 
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results.  The  methanol  absorptivities  for  reference  17  are  within 
9%  of  the  EVC  findings  for  the  P  and  R  bands.  The  methanol  P  and 
R  band  absorptivity  results  for  reference  14  differ  about  20% 
from  those  of  the  EVC.  The  remaining  absorptivities  in  Table  6 
vary  from  the  EVC  evaluations  by  30  to  70%.  This  clearly 
indicates  a  need  for  accurate  and  critically  evaluated  vapor 
phase  absorbance  measurements. 

The  EVC  method  provides  a  solute  vapor  pressure 
accurate  to  within  2%.  The  EVC  absorbance  versus  vapor  pressure 
measurements  provide  slopes  with  a  precision  of  ±2%  or  better. 
These  slopes  are  measured  over  the  Beer's  Law  range  for  a 
statistically  significant  number  of  concentrations.  In  the  case 
of  acetone  only  four  vapor  concentrations  show  adherence  to 
Beer's  Law  for  the  intense  1217  and  1365  cm'’  bands.  This  is 
shown  with  the  absorbance  versus  vapor  pressure  plot  in  Figure  5 
for  the  bands  of  Table  6.  An  added  benefit  of  the  EVC  method  is 
the  controlled  presence  of  water  vapor.  This  permits  the 
assessment  of  the  absorptivities  for  conditions  similar  to  those 
found  in  many  open-path  FTIR  spectrometry  applications. 


4 .  CONCLUSIONS 

Quantitative  measurement  of  vapors  in  the  atmosphere 
with  open-path  Fourier  transform  infrared  (FTIR)  spectrometry 
demands  the  use  of  accurate  absorptivities.  A  comparison  of 
interlaboratory  absorptivities  exhibits  significant  variation. 
The  integration  of  the  methanol  PQR  absorbance  band  furnishes  a 
criterion  for  correlating  absorbance  measurements  from  different 
laboratories. 

The  present  study  advocates  the  use  of  aqueous 
solutions  for  generation  of  selected  solute  vapors  at  specific 
concentrations  with  a  2%  accuracy.  Analysis  of  the  available 
literature  vapor  pressure  data  permits  an  assessment  of  the  best 
Wilson  coefficients  for  each  binary  system  studied.  The 
resultant  best  coefficient  values  along  with  the  vapor  pressures 
of  the  pure  solution  components  allow  an  accurate  evaluation  of 
the  solute  and  the  solvent  vapor  concentrations  for  equilibrium 
conditions.  An  air/vapor  mixture  is  continuously  closed  looped 
circulated  through  the  fixed  temperature  aqueous  solution  and 
infrared  (IR)  flow  cell  to  insure  solution  equilibrium.  A 
statistically  significant  number  of  vapor  concentrations 
establishes  the  linear  region  of  IR  absorbance  measurements 
(i.e.,  Beer  I s  Law).  The  resultant  experimental  precision  for 
the  absorptivities  is  ±2%  or  better  for  the  systems  studied. 
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Figure  1.  Equilibrium  Vapor  Cell 
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Table  2.  Antoine  Coefficients  for  Pure  Solution 
Components 


COEFFICIENTS 

COMPONENT 

A 

B 

C 

Water, 1 
Water, 2 

7.96681 

8.10765 

1668.210 

1750.286 

228.000 

235.000 

Acetone, 1 
Acetone, 2 

7.02447 

7.23967 

1162.000 

1279.870 

224.000 

237.500 

MEK(l) 

6.97421 

1209.600 

216.000 

Methanol , 1 
Methanol , 2 

7.87863 

8.07246 

1473.110 

1574.990 

230.000 

238.960 

Ethanol, 1 
Ethanol, 2 

8.04494 

8.16290 

1554.300 

1623.220 

222.650 

228.980 

Isopropanol,  1 
Isopropanol , 2 

6.66040 

7.75634 

813.055 

1356.142 

132.930 

197.970 
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Table  3.  Determining  of  Best  Wilson  Coefficients  for  Methanol-Water  System 
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Table  4.  Mean  and  Best  Wilson  Coefficient  Values 


Mean 

Best 

System 

w„ 

W21 

Vn 

Water-Methanol 

0.4271 

1.0833 

0.4085 

1.0900 

Water-Ethanol 

0.1395 

1.0260 

0.1950 

0.9104 

Water-Isopropanol 

0.0693 

0.7431 

0.0989 

0.6597 

Water-Acetone 

0.1814 

0.4896 

0.1768 

0.5148 

Water-MEK 

0.0835 

0.3415 

0.0242 

0.2721 

Table  5.  The  Effect  of  Apodization  on  Least  Squares  Slopes  for 
the  Methanol  Vapor  P  and  R  Branches  at  Various 
Resolutions.  All  slope  values  are  scaled  by  10^. 


{■ 

Boxcar 

APODIZATION 

Medium 

Norton-Beer 

Triangular 

1015.1 

4 

2.550±0.016 

2.551±0.015 

2.54310.016 

1014.8 

2 

2.712±0.021 

2.57810.016 

2.55210.016 

1013.7 

1 

3.260±0.025 

2.88710.018 

2.82610.018 

1013.7 

0.56 

3.200±0.019 

3.05410.018 

2.95010.017 

1057.6 

4 

3.751±0.029 

3.77710.029 

3.86110.037 

1057.3 

2 

3.742±0.027 

3.93010.029 

3.83610.037 

1055.2 

1 

3.914±0.021 

3.85110.050 

3.81410.023 

1057.6 

0.56 

3.963±0.033 

3.93910.032 

3.92010.032 

Table  6.  Comparison  of  Absorptivities  at  2  cm‘  Spectral  Resolution 
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